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ABSTRACT 

Understanding  turbulence  effects  on  wave  propagation  and  imaging  systems  has  been  an  active  research  area  for 
more  than  50  years.  Conventional  atmospheric  optics  methods  use  statistical  models  to  analyze  image  degradation 
effects  that  are  caused  by  turbulence.  In  this  paper,  we  intend  to  understand  atmospheric  turbulence  effects  using  a 
deterministic  signal  processing  and  imaging  theory  point  of  view  and  modeling.  The  model  simulates  the  formed 
imagery  by  a  lens  by  tracing  the  optical  rays  from  the  target  through  a  band  of  turbulence.  We  examine  the  nature 
of  the  turbulence-degraded  image,  and  identify  its  characteristics  as  the  parameters  of  the  band  of  turbulence,  e.g., 
its  width,  angle,  and  index  of  refraction,  are  varied.  Image  degradation  effects  due  to  turbulence,  such  as  image 
blurring  and  image  dancing,  are  revealed  by  this  signal  modeling.  We  show  that  in  fact  these  phenomena  can  be 
related  not  only  to  phase  errors  in  the  frequency  domain  of  the  image  but  also  a  2D  modulation  effect  in  the  image 
spectrum.  Results  with  simulated  and  realistic  data  are  provided. 

Keywords:  atmospheric  turbulence  effects,  turbulence-distorted  imagery,  signal  modeling,  Fourier  analysis, 
imaging  theory. 

1.  INTRODUCTION 

Atmosphere  turbulence  is  one  of  many  factors  that  affect  sensor  performance  and  image  quality.  It  results  in 
distorted  imagery  such  as  “image  blur”  and  “image  dancing”  degradations.  In  order  to  design  an  effective  signal 
processing  technique  to  remove  the  turbulence  degradation  phenomenon  or  improve  the  image  quality,  it  is 
important  to  understand  the  turbulence  effects  and  sources  of  degradation  first  from  a  signal  processing  and 
imaging  theory  point  of  view. 

Understanding  turbulence  effects  on  wave  propagation  and  imaging  systems  has  been  an  active  research  area  for 
more  than  50  years.  Conventional  atmospheric  optics  methods  use  statistical  models  to  analyze  image  degradation 
effects  that  are  caused  by  turbulence  [1-2].  These  methods  have  followed  the  theoretical  basis  from  contributions 
of  Kolmogorov  [3]  for  developing  a  statistical  model  for  the  spatial  structure  of  turbulent  air  flows,  Tatarskii  [4] 
for  applying  Kolmogorov’s  model  to  the  problems  of  wave  equation  and  propagation  through  regions  of  weak 
random  index  fluctuation,  and  Fried  [5]  for  extending  Tatarskii’s  results  to  describe  turbulence  effects  on  the 
resolution  of  imaging  systems. 

Image  degradation  effects  due  to  turbulence  include  image  blur  and  image  dancing  [1].  Conventionally,  many 
research  showed  that  turbulence-induced  wave  propagation  variations  across  the  aperture  of  an  imaging  system 
distort  the  point  spread  function  (PSF)  of  the  imaging  system.  The  influence  of  turbulence  blur  to  image  blur 
effects  has  been  characterized  in  terms  of  a  modulation  transfer  function  (MTF)  from  the  atmosphere,  that  was 
derived  by  Goodman  [6],  Hufnagel  et  al,  [7]  and  Fried  [8].  Later  work  on  the  model  for  turbulence  MTF  was 
extended  on  the  imaging  systems  in  the  visible  and  infrared  bands  [9-10].  The  evaluations  and  validations  of  the 
turbulence  MTF  in  the  infrared  imagery  were  conducted  in  Ref.  11-12.  However,  it  has  also  shown  [1,13]  that 
wave  propagation  through  turbulence  results  in  a  phase  error  in  the  frequency  domain  of  the  recorded  image.  One 
of  our  objectives  in  this  paper  is  to  characterize  the  image  blur  effect  in  terms  of  a  spectral  domain  phase 
degradation. 


Various  statistical  descriptions  of  the  optical  wave  through  random  media  are  also  used  to  characterize  the  image 
degradation  effects.  For  example,  the  angle-of-arrival  fluctuations  of  an  optical  wave  in  the  plane  of  the  receiver 
aperture  are  associated  with  image  dancing  in  the  focal  plane  of  an  imaging  system  [1].  The  beam  wander  was 
also  characterized  statistically  by  the  variance  of  the  beam  displacement  along  an  axis  or  by  the  variance  of  the 
magnitude  of  the  beam  displacement.  This  phenomenon  is  also  related  to  image  dancing. 

The  phase  screen  method  is  another  extending  development  to  model  the  propagation  of  an  optical  wave  as  a 
limiting  case  of  extended  turbulence  confined  to  a  thin  slab  screen  between  the  transmitter  and  receiver  [14]. 
There  are  a  number  of  approaches  to  generate  random  phase  screens  with  the  proper  point  statistics  and  spatial 
and  temporal  correlation  properties  [1-2,  14-15].  The  thin  phase  screen  models  permit  the  experimentalist  to 
predict  statistical  results  associated  with  laboratory  experiments.  However,  it  requires  critical  placement  of  the 
phase  screen  with  respect  to  transmitter  and  receiver  [1]. 

Intuitively,  the  image  dancing  degradation  is  caused  by  the  temporal  variations  of  the  turbulence-distorted 
imagery.  Unlike  conventional  statistical  models,  another  objective  in  this  paper  is  to  characterize  the  image 
dancing  effect  using  a  Fourier  spectral  analysis  of  recorded  images. 

In  this  paper,  we  intend  to  understand  atmospheric  turbulence  effects  using  deterministic  signal  modeling  and 
imaging  theory.  Our  model  simulates  the  formed  imagery  by  a  lens  by  tracing  the  optical  rays  from  the  target 
through  a  band  of  turbulence.  We  examine  the  nature  of  the  turbulence-degraded  image,  and  identify  its 
characteristics  as  the  parameters  of  the  band  of  turbulence,  e.g.,  its  width,  angle,  and  index  of  refraction,  are 
varied.  Image  degradation  effects  due  to  turbulence,  such  as  image  blur  and  image  dancing,  are  revealed  by  this 
signal  modeling.  We  show  that  in  fact  these  phenomena  can  be  related  not  only  to  phase  errors  in  the  frequency 
domain  of  the  image  but  also  a  2D  modulation  effect  in  the  image  spectrum. 

The  paper  is  organized  as  follows.  The  signal  propagation  model  and  imaging  without  and  with  turbulence,  for 
parallel  and  slanted  turbulence,  are  presented  in  Section  2.  This  signal  modeling  is  used  to  characterize  the  image 
blur  due  to  turbulence  in  terms  of  the  magnitude  and  phase  degradations  that  are  presented  in  Section  3.  Section  4 
presents  a  method  to  examine  the  spectral  properties  of  the  turbulence-distorted  images  that  are  related  to  image 
dancing  phenomenon.  The  results  of  magnitude  and  phase  degradations  of  simulated  data  are  presented  in  Section 
5.  Spectral  properties  of  real  field  data  are  presented  in  Section  6.  The  summary  is  provided  in  Section  7. 

2.  SIGNAL  MODELING 
2.1  Signal  Propagation  Model  and  Imaging 

Fig.  1  illustrates  the  geometry  of  image  formation  [16]  using  the  signal  propagation  modeling.  Assume  the  target 
is  a  planar  target  with  an  area  of  (-y  7 ).  Each  point  of  the  target  emits  a  time -dependent  signal  with  a  known 

frequency  co  or  wavelength  X  [  ®  =  (2x  c)  /  A  ],  c  is  the  speed  of  the  light,  that  is,  p(t )  =  exp(y®  t)  ■  We  examine  the 
signal  that  is  emitted  by  the  target  and  is  received  at  the  lens  with  a  size  of  u  e  (-Z,Z)  .  The  image  formed  by  the 
lens  is  calculated  using  a  wavefront  reconstruction  method  [17].  Asa  result,  the  image  distortion  is  due  to  the 
turbulence  only  since  the  wavefront  reconstruction  is  approximation  free. 


The  X  -coordinate  is  used  to  identify  range  from  the  target  to  the  lens,  and  V  specifies  the  cross-range  domain. 
The  lens  is  located  on  the  line  x  =  X  and  identified  by  (Xc,u) .  The  time  delay  of  the  signal  propagation  from 


(x,y)  in  the  target  region  to  an  aperture  point  on  the  lens  is 


c  c 


where  d(u)  =  Jxc2  +  («  -  y)2  is  the  distance  the  signal  travels  from  the  target  to  the  aperture  point.  The  total  signal 


received  at  the  lens  is  obtained  from  the  analysis  of  wavefront  reconstruction  imaging  method  in  [17]: 


(2) 


s(u,t)=  [  f(y,t)exp[jco t  +  j </>(}’, t )] exp[  y k d (w )] dy 

The  function  f(y,t )  is  the  reflectance  map  function  (related  to  the  desired  image)  that  contains  both  the  target’s 
physical  properties  and  fluctuations  of  the  amplitude  of  radiating  (emitting)  source;  and  ^(  yj)  is  an  unknown 
phase  function  that  represents  the  relative  phase  delay  of  the  imaging  wave  on  a  target  point  (related  to  the 
distance  between  the  source  and  the  point  target  and  the  unknown  phase  fluctuations  of  the  source),  and  the  lack 
of  coherence  among  the  target’s  radiating  points  and  the  receiver  structure  for  the  passive  array.  The  variable  k  is 
the  wavenumber,  k  =  co/c .  For  simplicity,  we  define 

g(y,  t )  =  /(.v,  t)  exp [j<f>(y,  0] 

Following  the  analysis  in  [17],  after  taking  the  spatial  Fourier  transform  of  both  side  of  (2)  with  respect  to  u ,  we 
obtain  Gv(ku,t)  =  exp (-jat ) exp [-j^k2  - k; Xc)S(ku,t)  (3) 

where  Gv(ku,t)  is  the  Fourier  transform  of  the  image  formed,  S(k,,j)  is  the  Fourier  transform  of  the  measured 
signal  at  the  lens,  and  ku  is  the  frequency  of  the  aperture  element  u  .  The  reconstructed  signal  g(y,t)  is  simply 
the  inverse  transform  of  Gr(ku,t )  • 

For  a  narrow  band  or  a  single  frequency  source,  the  reconstructed  signal  g(v,t)  can  be  denoted  as  the  function  of 
g(  v)  [15-16].  This  is  based  on  the  fact  that  the  recorded  image  is  the  integration  of  the  measurements  over  a 
period  of  time,  T  ,  e.g.,  the  camera  integration  time,  i.e., 

g{y)  =  \T\g{y,t)\dt 

For  a  wide  band  signal,  e.g.,  a  longwave  infrared  signal  with  a  bandwidth  of  8  - 12//  m  (micrometers),  the 
bandwidth  is  divided  into  a  number  of  narrow  wavelength  bands.  Then  the  signal  is  reconstructed  for  each 
wavelength  band. 

2.2  Imaging  in  Present  of  Turbulence 

Now,  we  consider  the  imaging  system  geometry  when  a  band  of  turbulence  of  width  D  is  present  in  Figs.  2-3. 

The  main  factor  that  has  to  be  incorporated  in  this  scenario  is  the  effect  of  turbulence  in  the  time  delay  of 
propagated  waves.  We  divide  the  distance  that  the  signal  travels  from  the  target  to  the  element  into  three  portions, 
d j ,  d2,  and  d3 ,  as  shown  in  Figs.  2-3.  In  this  case,  dx  represents  the  wave  propagation  path  from  the  target  to  the 

entry  point  of  the  turbulence,  d2  the  path  inside  the  turbulence,  and  ds  the  path  from  the  exit  point  of  the 

turbulence  to  the  lens.  The  total  time  delay  for  the  wave  propagation  is  obtained  as 

r  =  — +  — ,  v,=v3=c,  v2  =  c I n-i  (4) 

n  v2  v3 

where  n1  is  the  index  of  refraction  of  the  material  (substance)  in  the  turbulence  band. 

The  signal  received  at  the  lens  can  be  rewritten  as  follows: 

s(u,t)=  f(y,t)exp\jcot  +  j<f>(y,t)]exp\j(k1dl+k2d2+kld2)]dy  (5) 

where  k2  =kj  n  and  kx  is  the  wavenumber  of  the  air.  After  the  distances  of  dx,  d2 ,  and  d3  are  calculated,  the 
reconstructed  signal  is  obtained  the  same  way  as  in  (3). 

In  the  following,  we  determine  the  distances  of  dt ,  d^,  and  d3  using  ray  tracing.  We  consider  two  scenarios  of 

turbulence:  1)  parallel  turbulence  in  which  the  turbulence  is  parallel  with  the  target  and  sensor  (see  Fig.  2);  2) 
slanted  turbulence  in  which  the  turbulence  has  a  slanted  angle  with  respect  to  the  target  and  sensor  (see  Fig.  3). 


2.2.1  Parallel  Turbulence 


Fig.  2  shows  the  signal  propagation  geometry  when  a  band  of  turbulence  is  parallel  with  the  target  and  sensor 
plane.  The  signal  entrance  angles  in  the  air,  turbulence  band,  and  out  the  air  are  identified  as  Qx,  e2,  and  , 
respectively.  According  to  the  Snell  law,  we  have 

«,  sin0,  =  sin(?2  (6) 

where  nx  and  n2  are  the  index  of  refraction  of  the  air  and  the  turbulence,  respectively.  From  the  geometry  in  Fig. 
2,  we  have  following: 

yt  =  yi+yi+y^  (7) 

where  the  variables  yx,  y2,  y} ,  and  yt  are  the  distances  defined  in  Fig.  2.  According  to  the  signal  propagation 
geometry  shown  in  Appendix  A,  the  value  of  sin<9[  is  solved  by  minimizing  the  following  for  a  given  point  target 
at  ( x ,  v)  and  every  element  of  the  aperture  (that  is  u  ): 

nun(yl+y2+y3-yt) 

9, 

The  distances  of  dx,  d2,  and  d2  can  be  obtained  from  the  value  of  sin0, ,  the  target  range,  turbulence  width,  and 
the  distance  from  the  turbulence  to  the  lens  as  shown  in  Fig.  2. 

2.2.2  Slanted  Turbulence 

Fig.  3  shows  the  signal  propagation  geometry  when  a  band  of  turbulence  has  a  slanted  angle  0r  with  respect  to 
the  target  and  sensor  plane.  Now,  the  equation  of  (7)  is  still  satisfied  in  the  slanted  direction  with  angle  Q  as 
shown  in  Appendix  B.  Similarly,  the  value  of  sin0,is  solved  by  minimizing  the  following: 

mm(yl+y2+yi-y,) 

0t 

3.  MAGNITUDE  AND  PHASE  DEGRADATIONS 

Literatures  [1,  13]  show  that  turbulence  causes  a  phase  error  in  the  frequency  domain  of  the  recorded  image. 

Using  the  signal  modeling  described  in  the  previous  section  allows  us  to  study  the  nature  of  the  turbulence- 
degraded  images  in  terms  of  both  magnitude  and  phase  degradations.  The  magnitude  degradation  is  a  straight 
forward  analysis  by  the  above  mentioned  signal  modeling  and  is  demonstrated  in  Section  5.  The  analysis  of  the 
phase  degradation  due  to  turbulence  is  presented  in  this  section. 

Let  G(k  )  and  G,(Ay  )be  Fourier  transforms  of  g(y)  and  gt(  v) ,  respectively.  If  there  is  a  phase  difference 
<p(kr)  between  two  signals,  G(k ,)  and  G{(ky) ,  we  have: 

G(ky )  =  Gj  (ky )  exp  [j  <p(ky )]  (10) 

The  phase  difference  cp(k  )  can  be  obtained  from  calculating  the  phase  angle  of  the  correlation  function  Fc(kr  ) , 
that  is  obtained  as  follows: 

Fc(ky)  =  G(ky)G*l(k,)  =  \Gl(ky)\2  exp [-j(p{ky)]  (11) 

The  phase  angle  is  a  ID  function  with  respect  to  the  frequency  domain  ky  ■ 

4.  SPECTRAL  PROPERTIES 

Image  dancing  is  a  degradation  caused  by  the  temporal  variations  of  the  turbulence-distorted  imagery.  In  this 
section,  we  examine  the  spectral  properties  of  the  turbulence-distorted  images.  We  contemplate  that  what  is 


referred  to  as  “dancing”  or  “waving”  in  images  that  are  degraded  by  turbulence  is  in  fact  a  modulation 
phenomenon. 

As  it  is  known,  the  shift  of  the  spectrum  of  the  image  causes  the  modulation  in  the  spatial  domain.  The  IR  or 
visible  camera  is  a  non-coherent  imaging  system.  Thus,  the  recorded  image  is  a  full-wave  rectified  version  of  the 
reconstruction.  We  anticipate  to  see  the  properties  that  are  associated  with  the  spectrum  of  a  full-wave  rectified 
signal. 

To  study  the  spectral  properties  of  the  image,  we  calculate  the  3D  spectrum.  Consider  an  image  sequence 
be  f(x,y,t)  and  denote  its  3D  Fourier  transform  by  F(kx,kv,a>)  ■  The  ID  power  spectrum  with  respect  to  the  time 

domain  can  be  obtained  by  the  following: 

Si(®)  =  III  F(kx,ky,a>)  |2 

kx  ky 

Similarly,  the  power  spectrum  with  respect  to  the  X  domain  is  obtained  by  the  following: 

S2(*v)  =  II|mvA,,,®)|2 

k  co 

The  power  spectrum  with  respect  to  the  y  domain  is  obtained  by  the  following: 

kx  co 

We  also  consider  the  2D  power  spectrum  with  respect  to  the  time  domain  by  the  following: 

s4  (kx  ,ky)  =  ^  |  F(kx ,  ky ,  a>)f 

By  studying  the  power  spectrum  in  the  time,  X  ,  and  V  domains  separately,  we  can  examine  the  spectral 
properties  of  image  affected  by  turbulence.  This  helps  us  identify  the  problem  analytically  and  understand  the 
turbulence  phenomena. 

5.  MAGNITUDE  AND  PHASE  DEGRADATIONS  OF  SIMULATED  DATA 

In  this  section,  we  present  the  simulated  imagery  using  the  signal  modeling  method  described  in  Sections  2-3.  We 
examine  the  nature  of  the  turbulence-degraded  image  and  identify  its  characteristics  as  the  parameters  of  the  band 
of  turbulence  are  varied. 

There  are  three  types  of  the  parameters  in  the  signal  modeling  to  describe  the  lens,  the  target,  and  the  band  of 
turbulence,  respectively.  Referring  to  Figs.  2-3,  the  parameters  can  be  divided  into  three  groups  in  the  following. 

The  parameters  of  lens  include:  the  size  of  the  lens,  L  .  The  parameters  of  the  target  include:  the  location  of  the 
target,  (x,y) ;  the  area  of  the  target,  Y0 ;  and  the  wavelength  of  the  target,  X  .  The  parameters  of  the  turbulence 

include:  the  index  of  refraction,  n2 ;  the  width  of  the  turbulence,  D ;  the  slanted  angle,  0  ;  and  the  distance 
between  the  turbulence  to  lens,  xr0  .  Varying  above  parameters,  we  form  the  signals,  g(y)  and  g,(  v) ,  with  and 
without  turbulence,  respectively. 

5.1  Magnitude  Degradation  due  to  Turbulence 

5.1.1  Varying  Turbulence  Index  of  Refraction 

Fig.  4  shows  the  magnitudes  of  signals,  g(  v)  and  g,(  v) ,  with  and  without  turbulence,  respectively,  for  the  cases 
of  the  index  of  refraction  being  n2  =1.5  and  2. 


(12) 

(13) 

(14) 

(15) 


The  parameter  of  the  lens  is  L  =  0.2  m.  The  parameters  of  the  target  are  (x,  v)  =  (Xc,  0) ,  Xc  =  100  m,  Y0=  0.5  m, 
and  A  =  8  - 12//  m  (longwave  infrared).  We  divide  the  wavelengths  into  4  narrow  wavelength  bands.  Other 
parameters  of  the  band  of  the  turbulence  are  D  =  10  m,  Q  =  0°,  and  x  =  50  m. 

The  solid  lines  in  Fig.  4  represent  the  magnitude  of  signal  with  turbulence,  g(y) ,  the  dotted  lines  the  signal 
without  turbulence,  gj(y) ,  for  the  first  wavelength  band.  From  Fig.  4,  comparing  with  gj(  v) ,  g(y)  becomes 
more  smeared  when  the  index  of  refraction  increases.  This  indicates  that  larger  index  of  refraction  of  the 

turbulence  causes  more  degradation  of  the  signal. 

5.1.2  Varying  Turbulence  Width 

Fig.  5  shows  the  magnitudes  of  signals,  g(y)  and  g,(  v) ,  with  and  without  turbulence,  respectively,  for  the  cases 
of  the  width  of  the  turbulence  being  D  =  10  m  and  20m. 

The  parameter  of  the  lens  is  L  =  0.2  m.  The  parameters  of  the  target  are  (x,  v)  =  (X  ,  0) ,  Xc  =  100  m,  Y0  =  0.5  m, 
and  A  =  8  - 12//  m  (longwave  infrared).  The  wavelength  is  divided  into  4  wavelength  bands.  Other  parameters  of 
the  band  of  the  turbulence  are  n2  =  2,  9r  =  0° ,  and  x  „  =  50  m. 

The  solid  lines  in  Fig.  5  represent  the  magnitude  of  signal  with  turbulence,  g(v)  ,  the  dotted  lines  the  signal 
without  turbulence,  g,(y) .  From  Fig.  5,  comparing  with  gt(y) ,  g(y)  becomes  more  smeared  when  the  width  of 
the  turbulence  D  increases.  This  also  indicates  that  the  thicker  layer  of  the  turbulence  causes  more  degradation  of 
the  signal. 

5.1.3  Varying  Turbulence  Slanted  Angle 

Figs.  6-7  show  the  magnitudes  of  signals,  g(  v)  and  gj(y) ,  with  and  without  turbulence,  respectively,  for  the 
cases  of  the  slanted  angle  of  the  turbulence  being  0  =  10°  and  20° . 

The  parameter  of  the  lens  isZ  =  0.2  m.  The  parameters  of  the  target  are  (x,  v)  =  (X  ,  0) ,  X  =  100  m,  Y0  =  5  m,  and 
A  =  8  - 12// m  (longwave  infrared).  The  wavelength  is  divided  into  4  wavelength  bands.  Other  parameters  of  the 
band  of  the  turbulence  are  D  =  10  m,  «  =  2  ,  and  Xr0  =  50  m. 

The  solid  lines  in  Figs.  6-7  represent  the  magnitude  of  signal  with  turbulence,  g(  v) ,  the  dotted  lines  the  signal 
without  turbulence,  g,(  v) ,  for  the  first  wavelength  band.  From  Fig.  6,  the  signal  g(y)  has  a  shift  from  g,(y) . 
When  the  slanted  angle  of  the  turbulence  increases,  this  shift  increases  too.  Fig.  7  shows  the  signals  that  are 
zoomed  in  around  the  central  point.  Comparing  with  gj(v) ,  g(v)  appears  smeared.  This  indicates  that  the 
turbulence  not  only  causes  the  signal  blurred,  but  also  the  signal  shifted  when  the  turbulence  has  a  slanted  angle. 

If  one  is  targeting  an  object  when  the  turbulence  presents,  the  location  of  the  object  that  appears  in  the  image 
could  be  shifted  from  the  real  location  when  the  turbulence  appears  at  a  slanted  angle.  This  needs  to  be  corrected 
for  turbulence-distorted  imagery. 

5.2  Phase  Degradation  due  to  Turbulence 

In  this  section,  we  study  the  signal  phase  degradation  due  to  the  turbulence  using  the  method  presented  in  Section 
3.  The  parameter  of  the  lens  isZ  =  0.2  m.  The  parameters  of  the  target  are  (x,  y)  =  (X  ,  0) ,  X.  =  100  m,  Y0  =  0.5  m, 


and  X  =  8  - 12//  m  (longwave  infrared).  The  wavelength  is  divided  into  4  wavelength  bands.  The  fixed  parameters 
of  the  band  of  the  turbulence  are  0r  =  0°  and  v  =  50  m.  We  vary  the  turbulence  width  and  index  of  refraction  in 
the  following  2  cases:  1  )D  =  10m,  n2  =  1.5;  2)  D  =  20m,  n2=  2  ■ 

The  solid  lines  in  Fig.  8  show  the  phase  angle  curves,  <p(ky) ,  of  the  correlation  function  Fc(kv)  in  (11)  for  the 

above  2  cases  using  the  first  wavelength  band.  The  central  part  of  the  phase  angle  curve  appears  to  be  a  Gaussian 
function.  We  then  fit  the  Gaussian  function  to  the  phase  angle  curve.  The  Gaussian  function  of  the  phase  angle 
with  respect  to  the  frequency  domain,  k  ,  is  defined  as: 

kl 

<P(ky)  =  A  exp(-  \)  =  A  exp (-ak  ) 

Lo~ 

where  the  parameter  A  determines  the  amplitude,  a  =  l/(2cr) ,  which  is  inversely  proportional  to  the  standard 
deviation  of  the  Gaussian  function.  The  fitted  Gaussian  function  curves  are  plotted  as  dotted  lines  in  Fig.  8.  The 
parameters  (A, a)  that  are  used  to  fit  the  Gaussian  functions  for  the  above  2  cases  are  also  listed  in  Fig.  8. 

This  study  shows  that  the  phase  angle  due  to  the  turbulence  is  a  Gaussian  function  instead  of  a  quadratic  function. 
This  indicates  that  the  correction  of  the  image  blurring  caused  by  the  turbulence  due  to  the  phase  degradation 
could  be  realized  by  searching  the  parameter  space  of  (A, a)  ■ 

6.  SPECTRAL  PROPERTIES  OF  REAL  DATA 

A  sequence  of  60  frames  of  a  natural  scene  with  turbulence  is  selected.  Fig.  9(a)  shows  one  of  the  frames  from 
this  sequence.  The  spectrum  with  respect  to  the  time,  X ,  and  y  domains  are  shown  in  Figs.  9(b)-(d),  respectively. 
From  Fig.  9(b),  the  spectrum  shows  weak  harmonics.  That  indicates  that  a  minor  degree  of  waving  appears  in  the 
time  domain.  The  spectrum  with  respect  to  X  domain  that  is  shown  in  Fig.  9(c)  does  not  show  strong  harmonics, 
which  indicates  the  waving  does  not  appear,  or  not  easy  to  be  observed,  in  the  X  domain,  the  vertical  direction. 
This  is  indeed  the  case  from  viewing  the  video.  Fig.  9(d)  shows  the  spectrum  with  respect  to  y  domain  and  it 
shows  strong  harmonics.  This  indicates  the  strong  waving  in  the  y  domain,  the  horizontal  direction.  Indeed,  the 
video  of  the  sequence  shows  the  strong  waving  appearance  in  the  horizontal  direction. 

Figs.  9(e)-(f)  show  the  2D  power  spectral  properties  of  this  sequence.  The  2D  power  spectrum  with  respect  to  the 
time  domain  is  defined  in  (15).  The  distributions  of  the  power  spectrum  at  k  =0  and  k  =0  are  shown  in  Figs. 

y  x 

9(e)  and  (f),  respectively.  Fig.  9(e)  does  not  illustrate  prominent  harmonics.  Meanwhile,  Fig.  9(f)  shows  strong 
harmonics.  This  corresponds  to  the  fact  that  there  is  no  obvious  waving  appearance  in  the  vertical  direction  (  x ) 
but  waving  in  the  horizontal  direction  (  y  ). 

Next,  a  sequence  of  100  frames  of  a  test  bar  pattern  scene  with  turbulence  [18]  is  selected.  Fig.  10(a)  shows  one  of 
the  frames  from  the  sequence,  which  exhibits  a  strong  turbulence;  in  this  case,  the  test  bar  pattern  is  severely 
distorted.  The  spectrum  with  respect  to  the  time,  X  ,  and  y  domains  are  shown  in  Figs.  10(b)-(d),  respectively. 
Figs.  10(b)-(c)  do  not  show  prominent  harmonics.  The  video  of  the  sequence  verifies  that  there  is  no  obvious 
waving  in  the  time  domain  and  the  vertical  direction  ( X ).  Flowever,  strong  harmonics  can  be  seen  around  the 
middle  frequency  in  Fig.  10(d).  Indeed,  the  video  shows  a  strong  waving  in  the  horizontal  direction  (  y  ).  There 
are  also  some  fluctuations  near  the  zero  frequency  in  Fig.  10(d).  This  is  the  natural  look  of  the  spectrum  of  the  test 
pattern. 

Figs.  10(e)-(f)  show  the  2D  power  spectral  distribution  of  this  sequence.  The  distributions  of  the  power  spectrum 
at  kr  =  0  and  k ,  =  0  are  shown  in  Fig.  10(e)  and  (f),  respectively.  There  are  no  prominent  harmonics  illustrated 


in  Fig.  10(e).  Meanwhile,  strong  harmonics  are  shown  in  Fig.  10(f).  These  also  correspond  to  the  fact  that  there  is 
no  obvious  waving  in  the  vertical  direction  ( X )  but  there  is  strong  waving  in  the  horizontal  direction  (  y  ). 


Last,  a  sequence  of  100  frames  of  a  test  bar  pattern  scene  with  weak  turbulence  [18]  is  selected.  Fig.  1 1(a)  shows 
one  of  the  frames  of  this  sequence,  which  exhibits  a  weak  turbulence.  The  video  of  the  sequence  shows  no 
obvious  waving  in  the  time  domain  and  the  vertical  direction  ( X ),  and  weak  waving  in  the  horizontal  direction 
( y  ).  These  are  shown  in  the  spectrum  with  respect  to  the  time,  X ,  and  y  domains  in  Figs.  1 1  (b)-(d), 
respectively.  That  is,  there  are  no  obvious  harmonics  in  Figs.  1  l(b)-(c)  and  possible  weak  harmonics  in  the  lower 
frequency  band  in  Fig.  11(d).  In  Fig.  1 1(d),  harmonics  are  also  present  near  the  zero  frequency  of  the  spectrum. 
This  is  also  the  natural  look  of  the  spectrum  of  the  test  pattern  that  is  approximately  a  periodic  signal. 

Figs.  1 1  (e)-(f)  shows  the  2D  power  spectral  distribution  of  this  sequence.  The  distributions  of  the  power  spectrum 
at  k  =  0  and  kx  =  0  are  shown  in  Fig.  1 1(e)  and  (f),  respectively.  Fig.  1 1(e)  does  not  present  prominent 

harmonics.  Fig.  1 1(f)  shows  weak  harmonics.  This  also  reflects  the  fact  that  there  is  no  obvious  waving  in  the 
vertical  direction  (  X )  but  there  is  weak  waving  in  the  horizontal  direction  (  y  ). 

7.  SUMMARY 

In  this  paper,  we  presented  a  signal  modeling  method  to  study  the  nature  of  a  turbulence-degraded  image.  The 
proposed  signal  modeling  was  used  to  identify  a  turbulence-distorted  image’s  characteristics  as  the  parameters  of 
the  band  of  turbulence,  e.g.,  its  width,  angle,  and  index  of  refraction,  are  varied.  The  magnitude  and  phase 
degradations,  and  spectral  properties  of  the  simulated  data  and  field  trial  imagery  data  with  turbulence  were 
demonstrated.  These  include  image  degradations  due  to  turbulence,  such  as  image  blurring  and  image  dancing 
phenomena.  The  results  illustrated  that  these  phenomena  can  be  associated  not  only  to  phase  errors  in  the 
frequency  domain  of  the  image  but  also  a  2D  modulation  effect  in  the  image  spectrum.  Further  research  on 
temporal  spectrum  of  images  will  be  useful  to  understand  and  characterize  the  turbulence-distorted  imagery.  This 
study  provides  a  foundation  for  the  development  of  algorithms  to  correct  turbulence  distortion  in  images.  This  is 
achieved  by  identifying  an  analytical  model  for  turbulence  distortion  in  acquired  imagery.  This  study  also 
provides  a  tool  to  assess  the  severity  of  turbulence  in  imagery.  As  a  result,  the  user  has  a  tool  to  determine  the 
usefulness  of  the  measured  data. 
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Appendix  A.  Signal  Propagation  Geometry  with  a  Parallel  Turbulence 

In  Fig.  2,  x,  represents  the  distance  from  the  target  to  the  turbulence,  D  the  width  of  the  turbulence  band,  and  * 
the  distance  from  the  turbulence  to  the  array.  Rewrite  (7)  by  using  (6)  and  replacing  ,  y2 ,  and  y  with  X{ ,  D , 
x3 ,  0X ,  and  02 ,  we  have 

x\  iand\  +Dn\an6i  +X3  tan0j  —yt  =  0  ,  where  n  =  1  /  n2  ■ 

Using  the  single-variable  bounded  nonlinear  function  minimization  method,  the  value  of  sin6(  is  solved  by 
minimizing  (8),  that  is, 

min  (x\  tan  0\+Dn  tan  9\  +  X3  tan  0\-yt)’  w^ere  xi  =  x-D  -  xr0  and  x3  =  xr0 
8t 

Appendix  B.  Signal  Propagation  Geometry  with  a  Slanted  Turbulence 

In  Fig.  3,  0  represents  the  slanted  angle  of  the  turbulence  with  respect  to  the  target  and  sensor  plane.  The 
distances  of  yx ,  y2,  and  y3  are  calculated  as  follows: 

y,  =(x  cos  0r  +  y  sin  0r  -  xr0  -  D)  tan  0X 
y,  =  D  n  tan  0X 
y3  =  (x,.0  +  u  sin 0r) tan 0X 
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Figure  1.  Image  formation  geometry. 


Figure  2.  Image  formation  geometry  with  parallel  turbulence. 
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Figure  3.  Image  formation  geometry  with  slanted  turbulence. 


Figure  4.  Varying  the  turbulence  index  of  refraction.  The  solid  lines  represent  the  magnitude  of  signal 
with  turbulence,  g(  v)  ,  the  dotted  lines  the  signal  without  turbulence,  gt(y)  ■  This  indicates  that  larger 
index  of  refraction  of  the  turbulence  causes  more  degradation  of  the  signal. 


Figure  5.  Varying  the  turbulence  width.  The  solid  lines  represent  the  magnitude  of  signal  with  turbulence, 
g(  v)  ,  the  dotted  lines  the  signal  without  turbulence,  g,(v)  •  This  also  indicates  that  the  thicker  layer  of  the 
turbulence  causes  more  degradation  of  the  signal. 


Figure  6.  Varying  the  turbulence  slanted  angle.  The  solid  lines  represent  the  magnitude  of  signal  with 
turbulence,  g(  v) ,  the  dotted  lines  the  signal  without  turbulence,  gj(  v)  ■  A  turbulence  slanted  angle  causes 
a  shift  between  two  signals.  When  the  slanted  angle  of  the  turbulence  increases,  this  shift  increases  too. 
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Figure  7.  Varying  the  turbulence  slanted  angle  (con’t).  Zoomed  in  version  of  Figure  6.  The  solid  lines 
represent  the  magnitude  of  signal  with  turbulence,  g(y) ,  the  dotted  lines  the  signal  without  turbulence,  gj(y)  • 


Figure  8.  Phase  angle  of  the  correlation  function  for  varying  D  and  n2  •  The  solid  lines  represent  the  phase  angle  of 
the  correlation,  dotted  lines  the  fitted  Gaussian  function  curves,  and  the  parameters  (A,  a)  used  to  fit  the  Gaussian 
functions. 
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Figure  9.  Spectral  and  2D  power  spectral  properties  of  a  natural  scene. 
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